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ABSTRACT 



We present a statistical analysis of the final lens sample from the Sloan 
Digital Sky Survey Quasar Lens Search (SQLS). The number distribution of 
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a complete subsample of 19 lensed quasars selected from 50,836 source quasars 
is compared with theoretical expectations, with particular attention to the se- 
lection function. Assuming that the velocity function of galaxies does not 
evolve with redshift, the SQLS sample constrains the cosmological constant to 
^a = 0.79^o;Q7(stat.)^o;Qg(syst.) for a flat universe. The dark energy equation 
of state is found to be consistent with w = —1 when the SQLS is combined 
with constraints from baryon acoustic oscillation (BAO) measurements or re- 
sults from the Wilkinson Microwave Anisotropy Probe (WMAP). We also ob- 
tain simultaneous constraints on cosmological parameters and redshift evolu- 
tion of the galaxy velocity function, finding no evidence for redshift evolution at 
z < 1 in any combinations of constraints. For instance, number density evolu- 
tion quantified as v n = dln^/dh^l + z) and the velocity dispersion evolution 
u ff = dln<7*/(iln(l + z) are constrained to v n = l.OGJ^g^stat.^oJ^syst.) and 
u ff = — 0.05^g;}g(stat.)^Q;Q3(syst.) respectively when the SQLS result is combined 
with BAO and WMAP for flat models with a cosmological constant. We find 
that a significant amount of dark energy is preferred even after fully marginaliz- 
ing over the galaxy evolution parameters. Thus the statistics of lensed quasars 
robustly confirm the accelerated cosmic expansion. 



Subject headings: cosmological parameters - - cosmology: theory 
structure, evolution — gravitational lensing: strong 



galaxies: 



1. Introduction 



The statistics of gravitationally lensed quasars serve as a unique probe of cosmology 
(ITurner et al.lll984l ). An advantage of gravitational lensing is the simple and robust theoret- 
ical basis, which enables solid predictions of lensing signals for a given mass distribution. For 
instance, the probability that a distant quasar is strongly lensed by an intervening massive 
galaxy is known to be sensitive to the value of the cosmological constant or dark energy 
(JTurnerl Il990l : iFukugita et al.l Il990l ) . This result arises because the lensing probability is 
sensitive to the cosmological volume in the intermediate redshift range of z ~ 0.2 — 1 that is 
a sensitive function of dark energy. Aft e r many early studies based on sm a ll samples (e.g., 
Fukugitaet"al~1ll992l : iMaoz fc Rixlll993l : iKochaneki Eiilxi iFalco et al.l ll998J : Ichiba fc Yoshii 
19991 ) , recen t results are broadly consistent with the standard A-dominated flat co smological 
model (e.g. Jchae et alJliool JMitchell et al.ll2005l : Ichaell2007i : lOguri et al.ll2008ah . Unfortu- 
nately, small number statistics remain a significant limitation for the cosmological results. 



In addition to the number of lenses, the redshift distribution of lens galaxies contains 
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complementary information on cosmological parameters (jKochanekl 119921 ) . The differential 
probability distribution of lens redshifts is fairly insensitive to both the source quasar popu- 
lation and magnification bias. On the other hand, the lens redshift distribution does depend 
sensitively on any redshift evolution in the velocity function of galaxies, and ther efore it has 
been used to constrain both cosmological parameters an d galaxy evo l ution (e.g., lOfek et al. 



2003; Chae & Mao 2003; Matsumoto fc Futamase 2008; Chae 2010; Cao&Zhu 



2012 ). A 



caveat is that sample incompleteness, including bi ased sampling in the image separation and 
redshift spaces, can bias the results significantly ( jCapelo fc Natarajanl 120071 ). 



Statistical samples of strongly lensed quasars have been construc ted at both radio 
and optical wavele ngths. The Cosmic Lens All-Sky Survey (CLASS; iMyers et al.l 12003 



Browne et al.l 120031 ) represents the largest survey of strongly lensed quasars in the radio. 
CLASS identified 22 lens systems from a sample of ~ 16500 flat-spectrum radio sources, 
of which 13 lenses from ~ 9000 radio sources constitute a statistically well-defined sample. 
While the effects of dust extinction and foreground lens galaxies on the selection of strong 
lens systems are negligible in the radio, th e poorly characterized redshift distribut ion of the 
source population is a major problem (e.g.. lMunoz et al.ll2003t iMcKean et al.ll2007l ). Any un- 
certainties in the mean redshift of the source population directly translate into a systematic 
error on any dark energy constraint derived from the statistics of the CLASS lenses. 



Our lens survey, the Sloan Digital Sky Survey Quasar Lens Search (SOLS; lOguri et al. 



20061 . l2008al ; llnada et al.ll2008al |2010| . 120121 ) . is constructed at optical wavelengths and rep- 
resents the largest survey for gravitationally lensed quasars to date. The SQLS is built o n 
the spectroscopic quasar catalog of the Sloan Digital Sky Survey (SDSS; lYork et al.l 120001 ) . 
The redshifts are known for all the source quasars, allowing robust comparisons of the ob- 
served lensing probability with theoretical predictions to extract cosmological information. 
In addition, the spectroscopic data, and the rich color information of the 5-band SDSS imag- 
ing data, make the identification of the lens candidates quite efficient. While the relatively 
poor seeing of the SDSS images prevent us from discovering quasar lenses with small image 
separation (< 1"), we can make detailed simulations of lenses to characterize the selection 



function of our lens survey (lOguri et al.ll2006l . hereafter Paper I) 



In this paper, we pr esent a detailed statistical analysis of the final lens sample of the 
SQLS (llnada et al.ll2012l hereafter Paper V) based on SDSS Data Release 7 (DR7). The 
statistical sample consists of 26 quasar lenses selected from 50, 836 source quasars in the 
redshift range 0.6 < z < 2.2 with Galactic extinction corrected (ISchlegel et al.lll998l ) mag- 
nitudes brighter than i = 19.1. Paying particular attention to selection effects, we compare 
the number distribution of the strong lenses with model predictions to obtain cosmological 
and astrophysical information. Previous studies of strong lensing statistics have focused on 
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either constraining cosmological parameters with a fixed galaxy evolution model, or vice 
versa. In this paper we consider simultaneous constraints on cosmological parameters and 
galaxy evolution parameters with several different priors on the cosmological parameters 
from other cosmological probes. 

This paper is organized as follow. After briefly describing the lens sample used for 
the statistical analysis in £|2j we present our theoretical model in detail in §|3j Section [4] 
presents our main results on cosmological constraints. Constraints on redshift evolution of 
the velocity function are given in §0 We discuss the results in §0 and we conclude in §0 
We denote the present matter density as Qm, the dark energy density as ^de, the equation 
of state of dark energy as w (which is assumed to be a constant throughout the paper), and 
the normalized Hubble constant as h = if /(100kms _1 Mpc" 1 ). For the special case where 
dark energy is assumed to be a cosmological constant (w = —1), the dark energy density is 
denoted as ^a- 



2. Lensed Quasar Sample 



The SDSS (SDSS-I and SDSS-II Legacy Survey; lYork et al.l 120001 ) is an imaging and 
spectroscopic survey covering 10,000 sq uare degrees of the sky, using a dedicated wide-field 
2.5-meter telescope (JGunn et al.l 120061 ) at the Apac he Point Observatory in New Mexico 



USA. Images ta ken in five broad-band filters (ugriz, iFukugita et al.lll996l ; iGunn et al.lll998 



Doi et al.ll2010l ) are red uced with an auto mated pipeline, leading to an astrometric accuracy 
better than about O'.'l (jPier et al.l 120031 ). and a photometric zeropoint accuracy of about 



0.01 magnitude over the entire survey area (iHogg et al.ll200ll ; ISmith et al.ll2002l ; llvezic et al. 



2004 ; iTucker et al.ll2006t iPadmanabhan et al.ll2008l ). In addition to imaging, SDSS conducts 
spectroscopic observations with a multi-fiber spectrograph covering 3 800 A to 9200 A, wit h 
a resolution of R ~ 1800 for targets s elected by the imaging data dBlanton et al.l 12003). 



All the data are now publicly available flStoughton et al.ll2002t lAbazajian et al. 
20051 . 120091 : lAdelman-McCarthv et alJbood . I2OO7L l2008h . 



2003 



2004 



The final statistical lens sample of the SQLS is base d on the SDSS DR7 qu asar catalog 
containing 105,783 spectroscopically confirmed quasars ( jSchneider et al.l 120101 ). While the 
quasars are selected using several different selection criteria, a subsample of quasars with 
0.6 < z < 2.2 and Galactic extinction corrected point spread function (PSF) magnitudes of 
i < 19.1 is nearly complete, with the completeness indep endent of the SDSS image morphol- 
ogy (IRichards et al.ll2002L l2006t IVanden Berk et al.ll20051 ). The SQLS constructs a statistical 
lens sample based on this subsample of 50,836 quasars. We identify strong lens candidates 
from SDSS images of these sources using two selection algorithms. One examines the image 
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morphology to select smaller separation lenses (although still with 9 > 1"), and the other 
identifies companions of similar colors to select lens candidates at wider image separations 
out to 9 < 20". We also require that the magnitude difference between the lensed images is 
smaller than 1.25 mag, because lenses with larger magnitude differences are difficult to find 
in the SDSS images. Our simulations of SDSS images have shown that our selection algo- 
rithms are nearly complete over the chosen image separation (1" < 9 < 20") and magnitude 
difference (|Am| < 1.25) ranges (Paper I). The lens candidates are then observed with vari- 
ous facilities to determine if they are real lens systems or not. The final lens sample consists 
of 26 strongly-lensed quasars satisfying the criteria described above. Paper V supplies more 
details of the definition of the source quasar sample, the selection of lens candidates, and 
the subsequent observations. 

The 26 quasar lenses in the final SQLS statistical lens sample are summarized in Table [TJ 
Most of the lens redshifts are measured spectroscopically, either from direct measurement of 
the spectrum of the lens galaxy or using strong absorption lines in the quasar spectrum in 
which the redshift of the absorber matches the inferred lens redshift from the color and the 
apparent magnitude of the lens. For lenses without spectroscopic redshifts, we adopt the 
lens redshifts inferred from these colors and magnitudes and include their errors as described 
below. 

We estimate the /-band magnitude of the quasar components, Iqso, from the follow-up 
high angular-resolution imaging. We define Iqso m analogy to the PSF magnitude in the 
SDSS data, but without the contribution of the flux from the lens galaxy. We derive Iqso 
from the magnitude of the brightest quasar image Ibri and the total magnitude J to t in follow- 
up imaging. In Paper I it was shown, using simulations of SDSS images, that the effective 
magnification factor of lens systems relevant for the PSF magnitude is given by 

A* = filhat + (1 - AOjUbri, (1) 

where 

£ = ~[l + tanh (1.76 -1.780)], (2) 

the image separation 9 is in arcseconds and /i to t and // Dr i are the total magnification and the 
magnification of the brightest image, respectively. Given this result, we can estimate Iqso 
by 

J QSO = -2.5log{/H0-°- 4Z - + (1 -fi)10-° AI ^}. (3) 

Both /qso and the /-band magnitude of the lens galaxy, J ga i, are summarized in Table [TJ 



As in lOguri et al.l (j2008al . hereafter Paper III) and llnada et al.l (j2010l . hereafter Paper 



IV), we consider a subsample of the statistical lens sample for our final analysis. First, we 
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restrict the image separation range to # m in < 6 < 6 m3iX with 9 V 



1" and n 



4", since we 



are interested in galaxy-scale lensed q uasars whose lens potentials are well approximated by 



singular isothermal dis tributions (e.g. jRusin &: Kochanekll2005l ; iKoopmans et al.ll2006l 12009 



Sonnenfeld et al.ll2012[ ). Second, we require that the lens galaxy is fainter than the quasar 
images in the /-band, Iqso < Igai- This condition is necessary to ensure that emission from 
the lens galaxies does not affect our lens candidate selection and the quasar target selection 
in SDSS. We find that 19 of the 26 lenses pass these additional selection criteria. 



3. Theoretical Model 

We mostly follow the methodology described in Paper III to calculate the lensing proba- 
bilities for the SQLS lens sample, although we include a number of modifications and updates 
for making more accurate theoretical predictions. 



3.1. Lens Potential 



Various observations of galaxy-scale strong lenses have convincingly shown that the 
radial mass distribution of lens gal axies is, on average, wel l described by the isothermal 
distribution with p(r) oc r~ 2 (e.g., iRusin fc Kochanekl 120051 ; IKoopmans et al.l 120061 . 12009 



Sonnenfeld et al.ll2012l ). As in Paper III, we adopt the elliptical extension of the singular 
isothermal sphere, the singular isothermal ellipsoid (SIE). The convergence of the SIE is 

1 n 1 1/2 



k(x,v) 



e E x(e) 



1 



e) 2 x 2 + y 2 



(4) 



where #e is the Einstein radius and e is the ellipticity of the mass distribution. The Einstein 
radius is related to the galaxy velocity dispersion a v by 

2 A s 



e h 



4vr(^ 
c 



D n 



(5) 



where D\ s and D os are the angular diameter distances from lens to source and from observer 
to source, respectively. 

The parameter A(e) is the velocity dispersion normalization factor for non-spherical 
galaxies. Co mputing Me) requires assu m ptions abou t the three-dimensional shapes of lens 
galaxies (e.g. JKeeton &: Kochaneklll998l ). IChad (120031 ) computed A(e) for two extreme cases, 
galaxies having either oblate or prolate shapes. In our calculation, we use the following fitting 
formulae, 

Aou(e) ~ exp (0.108^ + 0.180e 2 + 0.797e 5 ) , (6) 



for the oblate case, and 



Ap ro (e) 



0.258e + 0.827e fc 



(7) 



for the prolate case (e < e max ~ 0.71158). In our fiducial model we assume an equal number 
of prolate and oblate galaxies to compute the dynamical normalization 

A(e) = /obiAobi(e) + (1 - /obi)A pro (e), (8) 

with / bi = 0.5 as a fiducial value. We assume that the distribution of the ellipticity is 
described by a Gaussian distribution with peak e and standard deviation a e , but truncated 
at e = and 0.9. Based on the axis ratio distributions of early -type galaxies in the SDSS 
Jchoi et all 120071 : IPadilla fc Straussl I2OO8I : iBernardi etaD boioh . we adopt e = 0.25 and 
a e = 0.2 as our fiducial parameters, which are slightly different from the values of e = 0.3 
and a e = 0.16 adopted in Paper III. 

In addition to the main lens galaxy, we include a contribution from line-of-sight density 
fluctuations to the lens potential in the form of a constant convergence and shear. While 
the effect of the external convergence and shear on the total lensing probability is small, 
the effect depends strongly on the image separation, such that the convergence and shear 
can have a large impact on lensing p robabilities at larger images separations of 9 > 3" 
(jOguri et al.l l2005bl ; iFaure et al.l |2009j) . In t his paper we employ th e probability distribu- 



tions of convergence and shear presented by iTakahashi et al.l (120 111 ), which have been de- 
rived from ray-tracing in high-resolution iV-body simulations. A caveat is that the matter 
fluctuations toward quasar lens systems may be biased compared to the fluctuation along 
the random line-of-sight directions, because the massive gala xies that are typical of strong 



lens systems are known to reside in dense environments (e.g., iTreu et al.ll2009l ; IFaure et al. 



201ll ). Such correlated matt er in the vicinity of the lens galaxy can contribute to external 



convergence and shear (e.g. , iKeeton et al. 



2004 ; iMomcheva et al.ll2006l : ISuvu et al. 



201 



1997 

F 



Fassnacht et al. 



Holder fc Schechterl 120031 : iDalal fc Watson 



201ll ). Thus in computing the 



probability distributions of the convergence n ext and shear 7 ex t coming from the line-of-sight 
matter fluctuations, we assume a source redshift of z s = 2, which is higher than the mean 
redshift of our lens sample but is still within its redshift range. 

Given the lens po tential, we c ompute the lensing cross section o"i ens numerically using 
the public code glafic (jOgurill2010l ). Specifically we randomly create many sources for each 
set of parameters (e.g., the ellipticity, external convergence and shear), and estimate <Ti cns as 

9(L/fM) 



v\> 



C11S.J 



du- 



(9) 



over the source plane positions u where multiple images are produced and the flux ratio of 
faint to bright images is larger than 10 -a5 for doubles (the flux ratio cutoff of the statistical 
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lens sample). The index i indicates the number of multiple images, with i = 2 for two- 
image lenses and i = 4 for four-image lenses. The parameter fi is the magnification factor 
for each source position, which is computed using Equation (TTJ. The magnification factor 
and the quasar luminosity function <&(L) (see §3.3p are needed to include the effect of the 



magnification bias. The cross sections are computed in units of the Einstein radius #e and 
as a function of dimensionless image separation defined by 9 - 



3.2. Velocity Function 

The velocity function of galaxies is an essential part of the theoretical prediction of the 
lensing probability. In Paper III, w e adopted the ve locity function of early-type galaxies 



obtained from the SDSS DR5 data (Choi et a" 



function for galaxies of all types ( iBernardi et al. 



2007). In this paper, we use the velocity 
20101 ). rather than that of early-type galaxies 



only. One of the reasons for using the all-type velocity functions for the analysis lies in the 
difficulty in making robust morphological classifications of the lens galaxy population. We 
note, however, that our sample of strong lenses is dominated by early-type galaxies because 
of the image separation cut of 9 > 1". 

Our fiducial velocity function is derived from galaxies of all types with the velocity 
dispersion er„ > 125 km s" 1 in the SDSS DR6 data (IBernardi et al.ll2010f ). The functional 
form is 



dn 
da„ 



a v 

<7* 



exp 



a v 
a* 



P 



(3 da v 



T(a//3) a v 



(10) 



a 



0.41, and (5 = 2.59. For 



with = 2.611 x KT 2 (/i/0.7) 3 Mpc- 3 , a* = 159.6 km s" 
comparison, we also consider the velocity function of galaxies of all types presented by IChae 
( 20101). wh i ch is b ased on the velocity function measurements for the SDSS DR5 galaxies by 
Choi et all J2007h . 



dn 
da v 



Ov\ a T(a/(3) (a^ 
a J e T(a'/P) U* 



exp 



(3 da v 
T(a/P) a v 



111 



with = 7.4 x 10~ 2 /i 3 Mpc~ 3 , a, = 100.0km s" 1 , a = 0.69, = 2.01, a' = 6.61, and e = 0.044. 

These velocity functions were derived from the analysis of low-redshift z ~ 0.1 galaxies 
in the SDSS. We consider the possibility of redshift evolution of the velocity function by 
allowing 0* and cr* to evolve as power-laws with redshift, 



!• -)• (j)*{l + zf n , 



(12) 
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a*^a*(l + z) u % (13) 

where we approximated that values of 0* and a* for the two velocity functions given above 
are for z — 0. The case with v n = and u a = corresponds to the no evolution model that 
has been adopted in most of previous analyses of quasar lens statistics (e.g., Paper III). In 
what follows, we leave u n and v a as free parameters to be constrained by the data, except in 
84] where we assume no redshift evolution. 



3.3. Quasar Luminosity Function 

The quasar luminosity function (QLF) is needed to compute the magnification bias in 
Equation (jUJ). We adopt the latest quasar lumino sity function from th e combined analysis 



of the SDSS and 2dF data (2SLAQ) presented by ICroom et al.l (120091 ). Specifically we use 



the pure luminosity evolution model derived in the redshift range 0.4 < z < 2.3 of 

$(M g , Z) = lo0 . 4(1 _ /3h)(Ms _ M » ) - 1() o.4(l-A)(M 9 -Af*) ' ( M ) 

M*Az) = M*(0) - 2.5(kiz + k 2 z 2 ), (15) 



with parameters given by (/3 h , A, **, M*(0), ki, fc 2 ) = (3.33, 1.42, 1.45 xlO- 6 (/i/0.7) 3 Mpc mag 



9 



— 22.18 + 51og(/i/0.7), 1.44, —0.315). The luminosity function is given in terms of rest-frame 
g-band absolute magnitudes at z = 2 (i.e., M g = M g (z = 2)). We c onvert the QLF to ob 



served i-band apparent magnitudes using the K-correction derived in lRichards et al.l ( 120061 ). 
Given the fact that the quasar luminosity function was derived assuming Qm = 0.3 and 
^a = 0.7, we adopt these cosmological parameters for computing the absolute magnitudes 
used to compute the magnification bias in equation flU]) no matter what values of f2 M , ^de, 
and w we consider for the remainder of the analysis. 



3.4. Number of Lensed Quasars 

With the lensing cross section computed by Equation (jUJ), we compute the differential 
probability that a source at z = z s and with the SDSS i-band PSF magnitude i = «qso is 
strongly lensed by a lens galaxy at z = z\ with image separation 9 as 

d 2 pi . d 2 V f d6 da v dn 2 da\. 



d9dzi 



\ nta\ dV f d9 da v dn 2 da {cnSti 

Zs ^ so) = a{e) dz^dnJ Jd^d^ e *^r S{tQSo) > (16) 



where 9 E is the Einstein radius defined in Equation (jSJ), 6 = 9/9 E , Ci{9) is the completeness 
of our lens candidate selection estimated from simulations of the SDSS images (see Paper I, 
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note that Ci(6) ~ 1 for the range of 9 of our sample), 



d 2 V cdt 



dzidVt dz] 



[l + zi 



?D% 



and 



S(i 



QSOJ 



erfc 



^QSO — IgaiyVv^l) 



V2< 



<?i 



(17) 



'U 



The factor S(iqso) is inserted to take into account the condition that the lens galaxy must 
be fainter than the PSF magnitude of the quasar. The mean galaxy z-band magnitude i ga i is 
computed fr om the velocity dispe rsion <j v given the observed correlation (the Faber- Jackson 
relation) by iBernardi et al.l ( 120031 ). which in cludes the luminosity evolution of galaxies with 
redshift, together with K-correction from the lColeman et al.l (J1980J) template spectrum of an 
elliptical galaxy. While we consider all types of galaxies as lenses, we use this relation because 
early-type galaxies dominate the population of lens galaxies in our sample, especially because 
of the removal of small-separation lenses 9 < 1". A Heaviside step function was adopted 
for S(iqso) in Paper III, but here we add a Gaussian scatter. We use the typical observed 
scatter in the Faber- Jackson relation of <Tj = 0.5. 

Given the probability distribution in Equation ffTB]) . we can easily compute probability 
distributions as a function of the image separation or the lens redshift as 



dpi 

d9 

dpi 
dzi 



^s,tQSO) 



d d2pi 
i, d9dzi 



, &a i v 



s, ^QSOy 



d9 



d 2 Pi 
d9dzi 



and the total lensing probability is 



Pi{Z s , IQSO) 



d9 



dzi 



d 2 Pi 
d9dz] 



(19) 
(20) 

(21) 



The predicted total number of lensed quasars in our quasar sample is calculated by counting 
the number of quasars, weighted by the lensing probability. To speed up the computation, 
we follow the procedure introduced in Paper III to calculate the number of lensed quasars 
for each redshift-magnitude bin and then sum over the bins. We define the number of source 
quasars in the redshift range z s j — Az s /2 < z s < z s j + Az s /2 and a magnitude range 
^qso.fe — Az/2 < i qso < i qso ,k + Ai/2 to be N qso (z s j,i qso> i~). The number distributions and the 
total number of lensed quasars become 



d9 



Zs,j ^qso,fc 



1 'q80\%S,jl ^ 



qso,fcJ 



dpi 

~d~9 l 



z. 



s,ji ''qso.k) j 



(22) 
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dNi 

dzi 



/ v / j ^qso{Z s ,j 



;'; ^qso,fc, 



^s,j ^qso,fc 



dp i t 
dzi 



\%S,ji 1qso,k) j 



(23) 



and 



■iv i / J / u J- 'qso\Zs,ji ^qsOjk )Pi\Zs,ji 1qso,k)- v^^J 

^s,j ^qso,fc 

Again, the index i indicates the number of multiple images. The adopted bin widths of 
Az s = 0.1 and Az = 0.2 are same as those in Paper III. 



3.5. Likelihood 



We use the likelihood function introduced by iKochanekl (119931 ) to constrain model pa- 
rameters, 



ln£ 



Vi ( d2pi 

^ n \d6dz, 

lens 



(N 2 + N< 



4), 



(25) 



where d 2 pi/d9dzi is calculated from Equation ( fl6l) and N 2 (doubles) and N4 (quadruples) are 
from Equation (liH|) . An important improvement from Paper III is that we now include the 
lens redshift distribution as a constraint, particularly because lens redshifts are successfully 
measured for most of the lenses used for the analysis. However, there are three lens systems 
whose lens redshifts are still not determined well (see Table [p. For these lenses we include 
the lens redshift uncertainties assuming a Gaussian distribution, 



d 2 Pi 



[z s ,i QS o) -» / dzi 



exp 



*J 



zY 



2^ 2 



dQdzi " J ' V27TO-, 

where z and a z are listed in Table [TJ We use the estimator, 

Ax 2 = -21n(£/£ max ), 

to derive best-fit model parameters and confidence limits. 



d 2 Pi 
dOdzi 



y z s,lQSOj 



(26) 



(27) 



4. Constraints on Cosmological Parameters 

In this section, we constrain cosmological parameters by comparing the observed lensed 
quasars in the SQLS DR7 statistical sample with theoretical model predictions. For now 
we assume that the velocity function of galaxies does not evolve with redshift (i.e., v n = 
and v a = in Equations f|T2|) and ffT3|) ). although redshift evolution is considered when we 
estimate systematic errors. When necessary, the constraints are combined with those from 
baryon acoustic oscillation and cosmic microwave background anisotropy measurements. 
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4.1. Flat Models with a Cosmological Constant 

We start with the simplest model, a flat universe where dark energy is described as 
a cosmological constant Q\ (i.e., w = —1). This model has only one free parameter, Q\. 
Before deriving constraints on Q^, we compare the number distribution of lenses in our 
sample as a function of image separation with the model predictions, as shown in Figure [TJ 
The total number of lenses is indeed sensitive to Q\, and models with Q\ ~ 0.8 are broadly 
consistent with the observed number distribution. In Figure |2j we examine the normalized 
lens redshift distribution N~ l dN/dzi. Here we see that the observed lens redshift distribution 
is more consistent with models with smaller Q^ because of the relatively small number of 
lens galaxies at high lens redshifts. 

We compute the likelihood of Equation (1251) as a function of Q\. The result shown in 
Figure [3] indicates that the SQLS DR7 sample constrains the cosmological constant to be 
Q\ = 0.79^n'n7, where the error indicates the statistical la confidence limit. This is consistent 



'a — u - ' y_o.07' 



with our previous results presented in Paper III and llnada et al.l (120101 . Paper IV). The model 
with ^a = is rejected at 6a level. Dark energy is detected at high significance by the quasar 
lens statistics. 

If we divide the likelihood into the part contributed by the numbers and separations 
as compared to the lens redshifts, we see that there is some tension. The numbers and 
separations, which we used in our previous studies, favor somewhat higher Q\ than the lens 
redshifts. In Figure [31 we show the resulting likelihoods as a function of Qa based only on the 
observed numbers and separations of lenses, as well as only on the lens redshift distribution. 
While the tension is only at about the 2a level, this result may be suggestive of redshift 
evolution in the velocity function. This is one of the reasons that we consider simultaneous 
constraints on cosmological parameters and galaxy evolution in §[5j 



4.2. Non-Flat Models with a Cosmological Constant 

Next we relax the assumption of a flat universe, and consider cosmological constraints in 
the Qm-^a plane. Figure [H sho ws the constraint from the SQLS D R7 using the full likelihood 



model. As already known (e.g.. lKochaneklll996l ; IChae et al.ll2002l ). the degeneracy direction 



of lens constraints in the Qm-^a plane resembles that from Type la supernovae (see, e.g., 



Suzuki et al.ll2012l . for a recent result). We find that a cosmological constant is required at 



the 4<r level even for this non-flat case. 

Because of different degeneracy directions, cosmological parameters are better con- 
strained by combining several different cosmological probes. In this paper, we combine 
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our constraints with either those from the baryon acoustic oscillation (BAO) measurement 
or from the cosmic microwave background (CMB). The former uses the baryon wiggle 
in the matter power spectrum as a sta ndard ruler. We adopt results from the WiggleZ 



Dark Energy Survey (IBlake et all 120111 ). which measures the BAO scale at z — 0.6, com- 



bined with BAO me asurements in the SDSS luminous red galaxies at z = 0.2 and 0.35 



(jPercival et al.ll2010i ). We consider the anisotropy measured by the Wilkinson Microwave 



Anisotro py Probe (WMAP) a s the latter. Specifically we consider the seven-year WMAP 



result by lKomatsu et al.l (1201 ll ). and compute the likelihood for each cosmological parameter 
set by the so-called "distance prior" which encapsulates all the key information relevant for 
dark energy studies derived from the WMAP data. While our lensing constraints from the 
SQLS DR7 do not have any dependence on the Hubble constant, both the BAO and WMAP 
constraints are sensitive to the adopted Hubble constant. Thus, when adding the BAO and 
WMAP we always include the Hubble constant as an additional free parameter over which 
we marginalize to obtain constraints on parameters of interest. 

These BAO and WMAP constraints are also shown in Figure HI The best-fit parameters 
and la statistical errors are Qm = 0.281q;q3 and Q\ = 0.88to!io when the SQLS is combined 
with BAO, and Q M = 0.20+°;°^ and Q A = 0.78+°;^ when combined with WMAP. The three 
constraints are complementary in the sense that their degeneracy directions are quite different 
from one another and the combined constraints give considerably stronger constraints in 
the Qm-Qa plane than any of the individual constraints. That our lensing constraints are 
consistent with both the BAO and WMAP constraints is an important cross check of the 
current standard cosmological model. 



4.3. Flat Dark Energy Models 

Next we consider flat models where the dark energy equation of state w is allowed to 
vary. This model is parametrized by Qm(— 1 — ^de) and w. Figure [5] shows constraints 
from the SQLS DR7 as well as those from BAO and WMAP. The degeneracy direction of our 
lensing constraint in this plane is again similar to that of Type la supernovae, and hence is 
complementary to the BAO and WMAP constraints. The combination of these constraints 
suggests a cosmological model with Qm ~ 0.3 and w ~ — 1. 

We also show how Qm and w are constrained when the lensing information is combined 
with either the BAO or WMAP result in Figure [6] and Figured respectively. For the former 
case, the resulting constraints are Qm = 0.25ig!o3 an d w — ~ l-44lto*2§, an d f° r the latter 
case Qm = 0.23+g;o3 and w = — 1.19i^\ 7 . In both cases, Qm and w are reasonably well 
constrained thanks to the different degeneracy directions of the tests. 
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4.4. Systematic Errors 

Thus far we have considered only statistical errors. Our model involves several uncer- 
tainties and assumptions which act as systematic errors in our cosmological analysis. Here 
we estimate these systematic errors in a similar way as done in Paper III. Specifically we 
consider the following sources of systematic errors. 

• We computed the dynamical normalization A(e) in Equation (jlj) assuming galaxies 
consist of oblate and prolate populations with equal fractions / bi = 0.5. As in Paper 
III, we change the fraction by ±0.25 to estimate the uncertainty associated with this 
assumption. 

• The ellipticity distribution of lens galaxies is assumed to be a truncated Gaussian with 
a peak e = 0.25 and a width <j e = 0.2. We shift the peak e by ±0.1 without changing 
the dispersion in order to see how the ellipticity distribution affects the cosmological 
results. 

• We included the line-of-sight convergenc e K fiXt , and shear y e ^ t u sing the PDFs derived 
from ray-tracing of iV-body simulations ( jTakahashi et al.l 120 111 ) for a fixed source red- 
shift of z s = 2. We change the assumed redshift by ±1 to estimate the systematic 
error. 

• The faint end slope of the quasar luminosity function in Equation ( 1141) is needed for 
computing the magnification bias, yet current measurements of the slope are fairly 
uncer tain. Considering r esults of measurements on the quasar luminosity function 



(e.g., [Hopkins et al.l 120071 ). we change the faint end slope (3\ by ±0.1 while fixing the 
other parameters of the quasar luminosity functions. While the range of /3j considered 
here is smaller than what was adopt ed in Paper HI, it is still much larger than the 
measurement uncertainty reported in ICroom et al.l (120091 ). 



The velocity function given by Equation ([TO]) is another important source of system- 
atic error. W h ile w e adopted the velocity function measurement in the SDSS by 



Bernardi et al.l ( 120101 ) as our fiducial model, we investigate how t he cosmologi cal re- 
sults are altered by adopting the velocity function measurement by I Chad (J2010I ). The 
specific forms of both the velocity functions were given in §3.21 

While we made the assumption that the velocity function does not evolve with redshift, 
i.e., v n = and v a = in Equations (TT21 and (1T31) . respectively, we check the effect of 
redshift evolution by adopting the evolution of v n = —0.229 and v a = —0.01 predicted 
by the semi-analytic model of iKang et al.l ( 120051 ). Note that we will explore the effect 
of letting v n and v a be free parameters in § [51 
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• As discussed in Paper III, the condition zqso — igai < is arbitrary. To estimate the 
systematic error, we shift the condition to ^qso — z g ai < ±0.25, within which the lens 
sample used for the statistics does not change. 

In Table |2j we show the contribution of each of these systematic errors to the systematic 
error on Q\ for the flat models with a cosmological constant. We find that the largest sources 
of systemic error are the dynamical normalization / bi and the velocity function dn/da v , 
followed by the faint end slope of the quasar luminosity function. The finding is consistent 
with the earlier results in Paper III. The resulting final constraint on Q\ from the SQLS 
alone including the systematic error is Cl\ = 0.79io!o7( s ^ a ^-)-ao6( s y s ^-)- Table |3] summarizes 
the cosmological constraints for the three cases studied above, including our estimates of 
systematic errors. All the results remain consistent with the current standard cosmological 
model (Qm ~ 0.3, f2 DE ~ 0.7, and w ~ —1). 



5. Evolution of the Velocity Function 

Thus far we have concentrated on the constraints on cosmological parameters from 
the statistics of lensed quasars. However, the lensing statistics also allow us to study the 
evolution and structure of the galaxies that act as lenses. In particular, the lens statistics 
serve as a useful probe of the velocity function of massive galaxies at intermediate redshifts, 
0-2 < z < 1, which are difficult to observe directly. Indeed, the analysis in §4.41 suggests 
that redshift evolution of the velocity function is one of the most significant sources of 
systematic error. Here we relax the assumption about the absence of redshift evolution in the 
velocity function, and investigate redshift evolution parametrized by v n (Equation (TT21 ) and 
v a (Equation (TT31) ). Unlike previous studies, however, we still allow cosmological parameters 
to vary and consider simultaneous constraints on the galaxy evolution and cosmological 
parameters with the help of external cosmological probes such as BAO and WMAP. 

First we consider the flat models with a cosmological constant. With the two additional 
parameters v n and u a , this model now has three parameters. The results using SQLS alone 
are shown in Figure [HJ We find that a cosmological constant is still required even after the 
evolution parameters are left as free parameters and are fully marginalized over. Specifi- 
cally, the cosmological constant is constrained to be Ha = 0. 71^25 (stat.) l !io ( s y s ^0 from 
the SQLS alone, which is consistent with the result assuming no redshift evolution of the 
velocity function shown in §4.11 The model without cosmological constant (Q& — 0) is still 
inconsistent with the data at more than 2a. The constraint projected in the v n -v a plane also 
indicates that the SQLS data are consistent with the no-redshift-evolution case (u n — v a — 0) 
at IcT. Also note that the most degenerate direction in the evolution parameters roughly 
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corresponds to no evolution in the lensing optical depth. The additional constraints from 
BAO and WMAP, which significantly improve the constraint on Q^, similarly improve the 
constraints in the v n -v a plane. With these additional constraints, results are still consistent 
with the no-redshift-evolution case. Specifically, the measured values of the two parame- 
ters are v n = l.Q6±{H(st8Lt.)toil(syst.) and v a = -0.05+o;ig(stat.)+o;o3(syst.) when all three 
constraints are combined. The systematic errors are estimated in the same way as in §4.4[ 
including the sources of the systematic error other than redshift evolution of the velocity 
function. 

We conduct similar analyses of the simultaneous constraints in the non-flat models with a 
cosmological constant and in the flat dark energy models. The results, shown in Figures[9]and 
[TUl indicate that the cosmological constraints from SQLS become much weaker when we allow 
the evolution parameters to vary, although some useful constraints are still obtained. For 
instance, in the non-flat cosmological constant case, non-zero Q^ is again preferred at about 
2a. On the other hand, the evolution parameters are constrained reasonably well even after 
marginalizing over cosmological parameters. In all the cases, the SQLS data are consistent 
with no redshift evolution, which supports early (z > 1) formation and passive evolution of 
massive early-type galaxies. We give a summary of the constraints in Table HI Note that 
systematic errors on the cosmological parameters appear smaller than in Table [31 because 
of weaker cosmological constraints from the SQLS after marginalizing over the evolution 
parameters. 

Our non-detection of redshift e volution of t he ve l ocity function i s con si stent with earlier 



result s bv lMao fc Kochanekl ( 



19941 . lOfek et all (120031 ) . IChae & Mad (120031 ) . ICapelo & Nataraian 



(120071 ). and iMatsumoto fc Futamasd (120081 ). although we believe our result is more robust 
given the carefully controlled lens sample from the SQLS and the exp loration of the de- 
generacy with cosmological parameters. On the other hand, I Chad (120101 ) reported that the 
lens data imply redshift evolution of the velocity function once a more complicated model 
is adopted for the evolution. To check this result, we consider the evolution of the shape of 
the velocity function by replacing a and f3 in Equation (fTUl) as I Chad (120101 ) has suggested, 



a — > a I 1 + k 



e /8 



l + z 



m 



(3^(3[l + k 



c /3 



l + z 



(29) 



These additional parameterizations lead to differential redshift evolution in the number den- 
sity of galaxies with different velocity dispersions. In particular, redshift evolution naively 
expected from the redshift dependence of the h alo mass function, which predicts stronger red- 
shift evolution for larger velocity dispersions (IMitchell et al.ll2005l ; IMatsumoto fc Futamase 
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20081 ). can well be described by this parametrization (IChadl2010l ). 



Figure [TT] shows constraints in the v n -v a plane in the flat models with a cosmological 
constant using all three cosmological constraints and marginalizing over the additional evo- 
lution parameter kp as well as the cosmological constant Q\. We find no redshift evolution 
even in this case. The constraints on the individual evolution parameters are v n = 2A7~^l'jl, 
v a = —0.35^49, and kp = —0.35^35 (statistical errors only), which are fully consistent with 
the fiducial no-redshift-evolution case of u n = v ff = kp = 0. Our different conclusion might 
be due to the fact that our source quasar sample is restricted to z s < 2.2 and therefo re probes 
relatively low lens redshifts of zi < 1, whereas the lens sample used by I Chad ( 120101 ) extends 
to lens redshifts of zi > 1. On the other hand, our SQLS sample is more complete and has a 
better understood selection function th an the some what heterogeneous, incomplete sample 
of source and lens redshifts adopted by I Chad (120 10l ). 



Discussions 



6.1. Comparisons with Other Studies on Galaxy Evolution 



While our results are consistent with no evolution, the constraints are relatively weak. 
In the v n -v a plane, the evolution is roughly constrained to keep the lensing optical depth 
(oc 0*0"^) constant while weakly constraining the orthogonal combination. Studies of galaxy 
evolution have typically found significant number and mass evolution in the early-type pop- 
ulation from z — to 1 with a decline in the abundance of early-t ype galaxies by roughly 
a factor of 2 by z = 1 (e.g., iFaber et al.l 120071 ; iBrown et al.l 120071 ) . which corresponds to 
v n ~ — 1. This change in number is consistent with our results, but it also requires an 
increase in the characteristic velocity dispersion of 20% (i.e., v a = 0.25). For the numbers to 
decline, the mean mass as indicated by the velocity dispersion has to increase. Eliminating 
this degeneracy requires samples of lenses large enough to cleanly measure the evolution of 
the average image separation with redshift (see §6.21) . 

Our results on the evolution in the velocity function can also be compared with recent 
evolution measurements in the stellar mass function. For instance, from the examination of 
galaxy populations at z ~ 1 — 2 it has been shown that the stellar mass function of galaxies 
indeed evolves from z = to 1 , though the mass dependence on the form of redshift evolution 



has y et to be fully clarified (e.g.. lllbert et aU2010l ; lMatsuoka fc Kawarall2010HBrammer et al. 



201 ll ). These studies showed that the number density of galaxies for a given stellar mass 
range can evolve by a factor of ~ 2 from z = to 1, which is in fact compatible with our 
results given the large errors on the evolution parameters for our lensing analysis. While the 
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direct evolution measurement of the velocity function is very challenging, iBezanson et al. 



( 120 111 ) has recently measured the evolution of the velocity function from z — to 1.5 by 
taking advantage of a scaling relation between velocity dispersion, stellar mass, and galaxy 
structural properties, and found weak redshift evolution up to z ~ 1, which is consistent 
with our result. 



6.2. Relation between Velocity Dispersion and Image Separation 

Another possible source of uncertainties inherent to our analysis is the relation between 
velocity dispersions and image separations. While this uncertainty is partly taken into 
account in our analysis via the systematic error from the dynamical normalization, the true 
uncertainties can potentially be larger given complexities such as velocity anisotropies and 
the detailed luminosity profiles of galaxies. However, with a large number of lenses we can 
in principle determine the relation from the data, because the ensemble average of image 
separations is given by (9) oc al with a proportiona lity factor that is almost independent of 



cosmological parameters (see iKochanek 



1993L Il996h . 



To explore the possibility of calibrating the relation between velocity dispersions and 
image separations from the lensing data, we consider a simple parametric model in which 
the Einstein radius given by Equation (^ is modified to be 

#e = (1 + 2 4vt(^) 2 ^, (30) 

\ c J D os 

where the parameter £ parametrizes the relation such that £ = if the velocity dispersion 
exactly matches the one used for our SIE models. We consider the flat models with a 
cosmological constant, and obtain simultaneous constraints on Q^ and £. The result, shown 
in Figure d2l indicates that constraints on £ from the SQLS are degenerate with Q\. The 
marginalized constraints on each parameter are £ = 0.08^q'oq and £7a = 0. 62^33 (statistical 
errors only). Thus the cosmological constraints become much weaker, although models with 
a significant cosmological constant are still preferred. If we combine all three cosmological 
constraints, the resulting constraint is £ = 0.051o'q3. The slightly positive value reflects the 
fact that lens statistics with £ = favor slightly larger Q\ than the best-fit values from BAO 
and WMAP, and that the observed mean image separation appears to be slightly higher 
than the models predict (see Figure [1]). Models with £ = are consistent with the data to 
better than 2a. 
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Conclusion 



We have conducted a statistical analysis of the final sample of strongly lensed quasars 
from the SQLS (Paper V). A subsample of 19 lenses selected from 50,836 quasars has been 
used to derive constraints on various cosmological parameters as well as the redshift evolution 
of the velocity function of galaxies. We have derived cosmological constraints assuming no 
redshift evolution of the velocity function. For the flat models with a cosmological constant, 
we have found VL\ = 0.79to;oy(stat.)^o;Qg(syst.) f rom the SQLS alone, which is in concord 
with other cosmological constraints. The model with Q\ = is rejected at 6<r (statistical 
errors only), which represents a significant detection of dark energy independent of Type la 
supernovae or other cosmological probes. The systematic error is comparable to the la 
statistical error, suggesting the importance of careful studies of the systematics for robust 
cosmological constraints from lens statistics. We have also derived constraints on non-flat 
models with a cosmological constant and flat dark energy models, by combining the SQLS 
results with independen t constraints from BA O and WMAP, and obtained results consistent 



with other studies (e.g.. lKomatsu et al.ll201ll ). These constraints are summarized in Table [31 



We have also derived simultaneous constraints on the cosmological parameters and 
redshift evolution of the velocity function. We parametrize redshift evolution by two pa- 
rameters v n (Equation ( 1T2j) ) and v g (Equation (lT3"j) ). The SQLS data still prefer a dark 
energy dominated universe even after marginalizing over the evolution parameters, with 

^a = 0.7llo'25(stat.)l(uo( s y s t-) f° r ^ e ^ a ^ m °dels with a cosmological constant. We have 
found no significant evidence for redshift evolution in the velocity function, for example 
v n = 1.06+i;|(stat.)+°;^(syst.) and v a = -0.05+^(stat.)+^(syst.), when the SQLS re- 
sults are combined with BAO and WMAP in the flat models with a cosmological constant. 
We summarize the simultaneous constraints in Table HI The results remain consistent with 
no redshift evolution even if we consider evolution in the shape of the velocity function. 
A cautionary note is that because of the relatively low-redshifts of our source quasars, the 
SQLS sample probes galaxy evolution only at z < 1. It is of great importance to extend 
the lens statistics like the SQLS to higher quasar redshifts in order to study the number 
evolution of massive galaxies further, as well as for better constraints on dark energy. Fu- 
ture wide-field surveys such as Pan-STARRS and Large Synoptic Survey Telescope will dis- 
cover thousands of lensed quasars efficiently by taking advantage of time-domain information 



(jOguri fc Marshall! 120101 ). which should be helpful for advancing such applications. 
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Table 1. DR7 Statistical Sample from Paper V 



Object 



N„ 



zf 



*PSF 



'QSO 



Rcf. 



SDSS J0246-0825 


2 


1.686 


0.723 


17.76 


1.09 


17.80 


20.78 


1,2,3 


SDSS J0746+4403 


2 


1.998 


0.513 


18.71 


1.08 


18.71 


19.62 


4,5 


SDSS J0806+2006 


2 


1.538 


0.573 


18.89 


1.49 


18.43 


20.16 


2,6,7 


SBS 0909+523 


2 


1.378 


0.830 


16.17 


1.11 


15.94 


18.81 


3,8,9 


SDSS J0924+0219 


4 


1.523 


0.394 


18.12 


1.81 


18.40 


19.36 


3,9,10,11 


FBQ0951+2635 


2 


1.246 


0.260 


17.24 


1.10 


16.54 


19.66 


2,3,12 


SDSS J1001+5027 


2 


1.841 


0.415 


17.32 


2.86 


17.32 


19.63 


13,14 


SDSS J1021+4913 


2 


1.720 


0.451 


18.99 


1.14 


18.85 


19.82 


14,15 


SDSS J1055+4628 


2 


1.249 


0.388 


18.76 


1.15 


18.86 


19.73 


5,14 


PG1115+080 


4 


1.735 


0.311 


15.97 


2.43 


16.40 


18.91 


3,16,17 


SDSS J1206+4332 


2 


1.789 


0.748 


18.46 


2.90 


18.05 


19.51 


13 


SDSS J1216+3529 


2 


2.013 


0.55 ±0.05 


19.08 


1.49 


18.30 


20.31 


18 


SDSS J1226-0006 


2 


1.126 


0.517 


18.23 


1.26 


18.67 


19.71 


3,19,20 


SDSS J1313+5151 


2 


1.877 


0.194 


17.70 


1.24 


17.10 


17.49 


21 


SDSS J1335+0118 


2 


1.571 


0.440 


17.54 


1.63 


17.17 


19.40 


3,20,22 


SDSS J1353+1138 


2 


1.624 


0.25 ±0.05 


16.47 


1.41 


16.36 


17.80 


6 


SDSS J1405+0959 


2 


1.810 


0.66 


19.05 


1.98 


18.70 


19.70 


23,24 


SDSS J1455+1447 


2 


1.424 


0.42 ±0.1 


18.22 


1.73 


18.21 


18.51 


5 


SDSS J1515+1511 


2 


2.054 


0.742 


18.05 


1.95 


17.53 


20.02 


23 


Q0957+561 


2 


1.413 


0.36 


16.68 


6.17 


16.69 


17.11 


3,25,26 


SDSS J1004+4112 


5 


1.740 


0.68 


18.82 


14.72 


18.44 


18.82 


3,27,28,29,30 


SDSS J1251+2935 


4 


0.802 


0.410 


18.86 


1.79 


19.32 


18.43 


31 


SDSS J1320+1644 


2 


1.502 


0.899 


18.88 


8.59 


18.36 


21.59 


32 


SDSS J1330+1810 


4 


1.393 


0.373 


18.35 


1.76 


18.34 


17.84 


33 


SDSS J1332+0347 


2 


1.438 


0.191 


17.89 


1.14 


18.57 


18.06 


34 


SDSS J1524+4409 


2 


1.210 


0.320 


18.76 


1.67 


19.49 


18.33 


18 



Note. — The subsample of 19 lenses above the horizontal solid line is used for the final statistical 
analysis. See text for details. 

a Source (quasar) redshifts from the SDSS data. 

b Redshifts of lens galaxies. Those with errors are lens redshifts inferred from the color and mag- 
nitude of the lens galaxy. 

c The PSF magnitude in the SDSS «-band magnitude corrected for Galactic extinction. 

d Maximum image separation in arcsec. 
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Johnson /-band (Vega) quasar magnitude without correcting for Galactic extinction (see text for 
details) . 

f Johnson /-band (Vega) magnitude of the lens galaxy without correcting for Galactic extinction. 
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Table 2. Breakdown of Systematic Errors for the Flat Models with a Cosmological 

Constant 



Source Error on Q\ 



/ ob i -)■ ±0.25 
e-»±0.1 

z s for K ext and 7 cxt ->• ±1 

A ->■ ±o.i 

Different dn/da v 
dn/da v evolution 
^qso — ^gai — > ±0.25 -o!on 



+0.037 
-0.032 

+0.000 
-0.006 

+0.003 
-0.000 

+0.019 
-0.018 

+0.000 
-0.046 

+0.033 
-0.000 



Total 



+0.057 
-0.060 



Note. - - The total errors are estimated 
from the quadrature sum of all the system- 
atic errors. 
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Table 3. Constraints on Cosmological Parameters 



Model 


Data 


fiM 


^DE 


w 


flat £l\ 


SQLS 


= i-n A 


n 7Q+O.O6+O.O6 
u - ' J -0.07-0.06 


= -1 


non-flat Q\ 
non-flat Q\ 
non-flat Q\ 


SQLS+BAO 

SQLS+WMAP 

SQLS+BAO+WMAP 


n 90+0.03+0.02 

u -^°-0.03-0.02 
n 9n +0.08+0.07 
u - zu -0.06-0.07 
n 9Q+O.O2+O.OO 
u - ZJ -0.02-0.00 


fl QO+0.09+0.07 
u - oo -0.10-0.09 
n 7O+0.05+0.05 
u - '°-0.06-0.05 
n 71 +0.02+0.00 
u - ' i -0.02-0.00 


I~ 1 



flat firjE 
flat fipjE 
flat firjE 


SQLS+BAO 

SQLS+WMAP 

SQLS+BAO+WMAP 


n 9 1- +0.03+0.03 
u - zd -0.03-0.02 
n 9O+0.04+0.03 
u -^°-0.03-0.03 

n 90+0.02+0.01 

u - zo -0.02-0.01 


= 1 
= 1 
= 1 


— Qm 

— Qm 

— si m 


1 a a +0.22+0. 17 
-■-•^^-0.25-0.18 
1 iQ+0.17+0.14 
1 - la -0.17-0.15 
1 1 1 +0.14+0.08 
i - ii -0.17-0.10 



Note. - - The parameter values are followed by their statistical and systematic un- 
certainties, respectively. See the last paragraph of §[T]for the definitions of cosmological 
parameters. 



Table 4. Constraints on Cosmological Parameters and the Redshift Evolution of the Velocity Function 



Model 



Data 



ft 



M 



n 



DE 



W 



flat Qa 
flat ftA 
flat f^A 
flat J1a 



SQLS = 1 - ft/ 

SQLS+BAO =1-0/ 

SQLS+WMAP = 1 - ft/ 

SQLS+BAO+WMAP = 1 - ft 



A 



0.71 
0.69 
0.73 
0.71 



20+0.12 
25-0.10 
03+0.00 
04-0.00 
02+0.00 
03-0.00 
02+0.00 
02-0.00 



1.071} 



1.15 



+i 
i 

1.00+} 
1.061} 



61+0.27 
64-0.82 
36+0.33 
43-0.72 
35+0.30 
.41-0.63 
36+0.33 
39-0.64 



-0.05j# 

-0.051°; 

-o.oelo 1 ; 

-0.051°; 



20+0.04 
20-0.05 
19+0.03 
16-0.03 
19+0.04 
16-0.03 
19+0.03 
16-0.03 



non-flat ftA 
non-flat ftA 
non-flat ftA 



SQLS+BAO 

SQLS+WMAP 

SQLS+BAO+WMAP 



0.25 



i-o. 
-o. 

0.321°; 
0.301°,; 



06+0.01 
05-0.01 
11+0.04 
09-0.13 
02+0.00 
02-0.00 



1.02; 



-0. 
-0. 

0.691° 
0.711°; 



20+0.06 
28-0.03 
07+0.10 
08-0.03 
02+0.00 
01-0.00 



0.15 



+1 
1 

1.291} 
1.061} 



61+0.20 
39-0.34 
24+0.38 
.62-0.66 
22+0.32 
.24-0.65 



-0.14; 



-0. 
-0. 

-0.061°, 
-0.051°; 



20+0.08 
23-0.02 
21+0.02 
16-0.06 
15+0.04 
16-0.03 



flat ft DE SQLS+BAO 

flat ft DE SQLS+WMAP 

flat ft DE SQLS+BAO+WMAP 



0.211°; 



0.17; 



-o 

-0 

0.29l° 



07+0.01 
05-0.01 
15+0.06 
06-0.05 
02+0.00 
03-0.01 



1-ft 
1-ft 
1-ft 



M 



M 



1 7n +0.45+0.05 
■l./U_ 4o_o.lO 

1 ci +0.69+0.34 
" ± - o± -0.60-0.32 



M 



-1.01 



hO. 14+0.03 
-0.18-0.03 



-0.48 

-0.51 
1.06 



71+0.33 
56-0.29 
38+0.88 
.18-0.80 
18+0.28 
39-0.66 



-0.04" 



-0. 

-0. 

-0.071°, 
-0.051°; 



19+0.04 
19-0.03 
17+0.04 
19-0.03 
16+0.03 
14-0.02 



CO 



Note. - The parameter values are followed by their statistical and systematic uncertainties, respectively. See the 
last paragraph of § CD for the definition of cosmological parameters. 
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T3 
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~\ — I — i — rrn — i — i — i — i — i — i — I — i — i — i — i — i — i — i — i — i — I — i — i — i — i — i — r 



DR7 statistical sample 
50,836 source quasars 



n A =0.95 
n A =0.75 




IT I '1- 



[arcsec] 



Fig. 1 . — The histogram shows the image separation distribution of the strong lenses in the 
statistical lens sample used for our cosmological analyses. The subsample contains 19 lenses 
selected out of 50,836 source quasars, as summarized in Tabled) Lines show the theoretical 
predictions for three different values of the cosmological constant fi\ assuming a flat universe 
and no evolution of the galaxy velocity function. The vertical dotted line shows the 9 min = 1" 
lower limit for the image separations in the statistical lens sample. 
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t 1 1 r 



n 1 1 r 



fi A =0.95 



n =0.75 




Fig. 2. — The histogram shows the normalized lens redshift distribution for our lens sample. 
Lines show theoretical predictions for three different cosmological models, as in Figure [TJ 
For lenses with errors on the lens redshift (see Tabled]), we adopt the best estimated values 
in constructing the histogram. 



-35- 




n A =i-n M 



Fig. 3. — Likelihood distributions as a function of the cosmo logical constant Q\ for flat 
universes. The dotted and dashed lines show the separate likelihood distributions for fitting 
the numbers and image separations alone (dotted) and lens redshifts alone (dashed). 
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a 



1 - 












0.5 



n 



1.5 



M 



Fig. 4. — Constraints on Qm and Q\ for the non-flat models with a cosmological con- 
stant. Contours show la and 2a confidence regions from the three different cosmological 
probes: S QLS strong lens stat i stics (this paper), baryon acoustic oscillation (BAO) mea- 
surements (IPercival et al.l I201CH : iBlake et al.l 1201 ll ) , and the CMB anisotropy from WMAP 
(JKomatsu et al.ll201ll ). The dotted line indicates a flat universe with Qm + ^a — 1- The 
upper left shaded region indicates models with no big bang. 
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0.6 



Fig. 5. — Constraints on Qm — 1 ~~ ^de and w for the flat dark energy models. As in 
Figure HI the la and la constraints from the three different cosmological probes are shown 
by contours. The horizontal dotted line indicates a cosmological constant with w = —1. 
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~~\ 1 r 




0.6 



^M ■*■ ^DE 



Fig. 6. — Similar to Figure but with constraints only from the SQLS DR7 and BAO. The 
innermost contours show the combined constraint. 
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0.6 



^M ■*■ ^DE 



Fig. 7. — Similar to Figure but with constraints only from the SQLS DR7 and WMAP. 
The innermost contours show the combined constraint. 
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Fig. 8. — Simultaneous constraints on cosmological parameters and redshift evolution of 
the velocity function of galaxies in the flat models with a cosmological constant. Upper: 
Likelihood distributions as a function of cosmological constant fl\. The solid line shows the 
likelihood from the SQLS alone with no redshift evolution as in Figure |3j The dashed line 
is the likelihood from the SQLS alone after marginalizing over the evolution parameters u n 
and v a . The dotted line is the likelihood distribution marginalizing over u n and v a when 
the SQLS result is combined with the BAO and WMAP constraints. Lower: Constraints 
on redshift evolution in the v n -v a plane after marginalizing over Q\. The outer contours 
are from the SQLS only, and the inner contours show the combined constraints from SQLS, 
BAO, and WMAP. Dotted lines in the lower panel indicate no redshift evolution (i/ n = 
and v a = 0). 



-41- 
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a 



0.5 



0.5 - 



^ b 
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■6 



■2 
1/ 
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Fig. 9. — Simultaneous constraints on cosmological parameters and redshift evolution of 
the velocity function of galaxies in the non-flat models with a cosmological constant. The 
upper panel shows constraints in the Qm-^a plane after marginalizing over the evolution 
parameters v n and v„. The dotted line is a flat universe with Qm + ^a = 1- Lower panel 
shows constraints in the v n -v a after marginalizing over the cosmological parameters. In both 
panels, from outer to inner contours, we show constraints from the SQLS alone, SQLS plus 
BAO, and the combination of all three probes (SQLS+BAO+WMAP), respectively. Dotted 
lines in the lower panel indicate no redshift evolution [y n = and v a — 0). 
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Fig. 10. — As in Figure |9j but for flat dark energy models with w as a free parameter. 
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Fig. 11. — Constraints on redshift evolution of the velocity function. Here we include 
the additional parameter kg wh i ch de scribes redshift evolution of the shape of the velocity 
function as proposed by IChad (120 ldl ). We consider the flat models with a cosmological 
constant. The SQLS constraint is combined with BAO and WMAP. We show constraints 
in the v n -v a after marginalizing over the other parameters Q\ and kp. Dotted lines in the 
lower panel indicate no redshift evolution [y n = and v a = 0). The marginalized constraint 
on kp is k = -0.35+°;^. 
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Fig. 12. — Constraints on the cosmo logical constant Q\ and the parameter £ in Equation (1301) 
that parametrizes the relation between velocity dispersions and image separations, for flat 
models with a cosmo logical constant. Outer and inner contours show constraints from the 
SQLS alone and from the combination of all three probes (SQLS+BAO+WMAP), respec- 
tively. The horizontal dotted line indicates the fiducial value £ = assumed in the paper. 



